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We extend the concept of classicality in quantum optics to spin states. We call a state "classical" 
if its density matrix can be decomposed as a weighted sum of angular momentum coherent states 
with positive weights. Classical spin states form a convex set C, which we fully characterize for a 
spin-1/2 and a spin-1. For arbitrary spin, we provide "non-classicality witnesses". For bipartite 
systems, C forms a subset of all separable states. A state of two spins-1/2 belongs to C if and only 
if it is separable, whereas for a spin-1/2 coupled to a spin-1, there are separable states which do 
not belong to C. We show that in general the question whether a state is in C can be answered by 
a linear programming algorithm. 
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I. INTRODUCTION 

The question of the classicality of quantum states has 
regained interest with the rise of quantum information 
theory Stronger-than-classical correlations between 
different systems are an important resource for quan- 
tum communication protocols, and the existence of large 
amounts of entanglement has been shown to be neces- 
sary for a quantum computational speed-up [2, 3] . How- 
ever, even for a single system the question of classical- 
ity is important. Historically the question goes back to 
two seminal papers in quantum optics by Sudarshan and 
Glauber [1, HI, who introduced the Glauber-Sudarshan 
P-representation for the states of a harmonic oscillator. 
This representation allows to decompose the density ma- 
trix in terms of coherent states of the harmonic oscilla- 
tor. For a single coherent state, the weight function of 
the P-representation (called P-function in the following 
for short) reduces to a delta function on the phase space 
point in which the coherent state is centered, and the dy- 
namics of the P-function is exactly the one of the clas- 
sical phase space distribution. It has therefore become 
customary in quantum optics to consider states with a 
positive P-function as classical. Several other criteria 
can be derived from this requirement. Using Bochner's 
theorem for the Fourier transform of a classical probabil- 
ity distribution jl], Richter and Vogel derived a hierarchy 
of observable criteria based on the characteristic function, 
which are both necessary and sufficient for classicality f^. 
This led to a recent demonstration of the negativity of 
the P-function in a quantum optical experiment |8] . Ko- 
rbicz et al. realized a connection of the positivity of the 
P-function to Hilbert's 17th problem of the decomposi- 
tion of a positive polynomial Since the P-function 
for a continuous variable system can be highly singular, a 
lot of attempts to define classicality have been based on 
other quasi-probability distributions (lo| as well, notably 
the Wigner function [ill, ■ 

These quasiprobability distributions for the harmonic 
oscillator have analogs for finite-dimensional angular 
momentum states [l^. The Wigner function for finite- 



dimensional systems has received a large amount of at- 
tention, ranginfffrom questions of its most appropriate 
definition [la . [3, [H, [H, [l3|, over classicality criteria 
[l3, to the importance of its negativity for quan- 
tum computational speed-up [IS] (see also for further 
references concerning the historical development of the 
Wigner function for finite-dimensional systems). Sur- 
prisingly, the P-function for finite-dimensional systems 
has been much less studied, in spite of its attractive 
mathematical properties. The P-function for a system 
with a finite-dimensional Hilbert space (i.e. formally a 
spin system) allows to decompose the densit y in atrix in 
terms of angular momentum coherent states [2l|. It can 
always be chosen to be a smooth function, expandable 
in a finite set of spherical harmonic functions [l3|. In 
contrast to the case of the harmonic oscillator, questions 
concerning the existence of the P-function (or its na- 
ture as a distribution or worse) do therefore not arise. 
This idyllic situation is somewhat perturbed, however, 
by the fact, already observed in [2l|, that for a spin sys- 
tem a large amount of freedom exists in the choice of the 
P-function, as it depends on two continuous variables 
on the Bloch sphere, whereas the density matrix for a 
system with c?-dimensional Hilbert space is specified by 
<P — 1 real independent entries. 

In this paper we show that the existence of a P- 
representation of the state of a spin system with a posi- 
tive P-function is a meaningful concept which allows to 
define the classicality of states of finite-dimensional sys- 
tems in a natural fashion, completely analogous to the 
classicality of the harmonic oscillator states of the elec- 
tromagnetic field. We shall call the corresponding states 
"P-representable", or P-rep for short. The set C of P- 
representable states form a convex domain in the space of 
density operators, containing the completely mixed state 
in its interior. We show that, surprisingly, all states of a 
single spin-1/2 are P-rep, and obtain an analytical cri- 
terion for P-representability in the case of a spin-1. For 
bipartite systems, the set of P-rep states is a subset of 
the set of separable states. For two spins-1 /2 the two sets 
coincide, whereas already for a spin-1/2 combined with 
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a spin-1, there are separable states which are not P-rep. 
We also show that the problem of deciding whether a 
given state is P-rep can be solved numerically by linear 
programming. 

In the following we will first motivate and define P- 
representability, then study simple cases of small spins, 
introduce a variational approach that gives rise to a lin- 
ear programming algorithm, and finally have a look at 
composite systems. We also develop some necessary con- 
ditions for P-representability based on measurable ob- 
servables, which may thus serve as "non-classicality wit- 
nesses", an extension of the by now well-known concept 
of entanglement witnesses j22| • 



II. DEFINITION OF P-REPRESENTABILITY 

A. Coherent states 

We first set some notations following the lines of 
Angular momentum coherent states are defined as eigen- 
states of and n.J with eigenvalues j{j + 1) and j, 
respectively, where n is a unit column vector which spec- 
ifies the quantization axis with polar angle and azimuth 
if, and J is the familiar angular momentum operator with 
components Jx , Jy and Jz ■ The transpose of the column 
vector n reads 

n {0, ipf = (sin cos ip, sin sin ip, cos 0) . 

An angular momentum coherent state can be expanded 
in terms of the states \jm) quantized on the z axis as 
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The coherent states form a complete, although not or- 
thogonal, basis set of normalized states within the space 
of the eigenfunctions of with given j, and 



2j + l 
An 



sin 0d0d(j)\0ip) {0ip\ = Uj+i, 



(1) 



where I2J-1-1 is the (2j -I- l)-dimensional identity matrix. 
We shall use the shorthand a = {0, (p) and denote da = 
sin 0d0d(t). The coherent state \0ip) associated with the 
vector n will be denoted In) or la). 



B. P representation 

The P-representation of a density operator p is an ex- 
pansion over the overcomplete basis of coherent states. 
This expansion reads 



p— daP {a) \a) {a\ . 



(2) 



where the P-function P (a) is real and normalized by the 
condition 



trp 



daP (a) = 1 



(3) 



If P (a) is non-negative then p is a classical mixture of 
pure coherent states with probability density P{a), and 
can therefore be considered as classical. In this case we 
shall say that p is P-representable, or "P-rep" for short. 

This definition has to be made more precise consider- 
ing that P(a) is not uniquely determined by the density 
operator. To show this, consider the multipole expansion 
of p, 

P ^ PkqTkq, Pkq = trpf^Q, (4) 

K={)Q=-K 

tkq = E (-ir"+^c;^'?,™jj"ii)(jm2i (5) 



where CfJ?,„, are the Clebsch-Gordan coefficients as 

I , J '"'IJ "''2 

[231 • Expanding the P-function as a sum of spherical 
harmonics. 
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one obtains a one-to-one relation between the coefficients 
of the two expansions for < K < 2j , 



Pkq 



= PkqV4tt 



y/T{2j-K + l)Ti2j + K + 2) 



(6) 



If _ftr > 2j the Euler Gamma functions in the denomina- 
tor become infinite; consequently regardless of Pkq the 
respective p^Q will be zero. It means that the choice of 
such Pkq is totally arbitrary. However, non-negativity 
of a P(a) for one choice of Pkq with K > 2j may be 
absent for another choice. Here is a simple example. 
Let the density operator be a projector on a coherent 
state, p = \ao) {ao\. An obvious P-function in this case 
is 5{a — ao); it can be considered non-negative since it 
can be approached by a sequence of non-negative func- 
tions, like Gaussians with decreasing width. An alter- 
native choice however would be to drop all non-physical 
terms in P with K > 2j, replacing the i5-function by a 
finite linear combination 



P(a) 
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which is not non-negative for all finite j (its tail away 
from the maximum at a = ao oscillates around zero). 

In view of the non-uniqueness of P(a) we reformulate 
the definition of P-representability demanding that the 
condition P > must be fulfilled at least for one partic- 
ular P(a). Under this definition the pure coherent state 
p = |ao) (aol will be P-rep, which is intuitively reason- 
able. We are thus led to the following definition: 
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Definition 1 A density matrix p is called P-rep if it can 
he written as a convex sum of coherent states, i.e. as in 
Eq. Q with a non-negative function P{a). 

We will now derive some simple consequences of this def- 
inition. 



C. Consequences 

Let V be the vector space of (2j + l) x (2j + l) hermitian 
matrices. The scalar product {X, Y) — trX'^Y defines an 
operator norm ||X|| — VtrX^X on V. We denote by 
Af the subset of non-negative density matrices, and by 
C the subset of P-rep states. The boundaries of these 
sets are respectively denoted dJ^ and dC. The following 
statements follow immediately from the above definition: 

1. The totally mixed state Pq = 2j+i'^^j+'^ P-rep, 
which is readily seen from Eq. ^ taking P (a) = 
l/47r. 

2. The set C of P-rep states is the convex hull of the 
set of coherent states. In particular, it is a convex 
set. 

3. Since all P-rep states are non-negative (but not 
vice versa) we have C C A/" C V. 

4. According to Caratheodory's theorem on convex 
sets applied to the {2j + l)^-dimensional vector 
space V, any non-negative Hermitian matrix can be 
represented as a convex sum of at most (2j -1-1)^-1-1 
projectors onto coherent states. In the case of 
density matrices subject to the condition trp = 1 
this number is decreased by 1. Finding a P- 
representation for a state p is thus equivalent to 
finding real non-negative coefficients Xi and coher- 
ent states \ai) such that 

(2j + l)' 

1=1 

5. A pure state is P-rep if and only if it is a coherent 
state. 

Proof The "if" part is trivial. For the "only if" 
part, assume that a state p is P-rep, i.e. that there 
exists a decomposition such as in We have 

= < (Ei^O^ 1, where 

equality occurs only for Kaijaj)! — 1 for all 
The latter condition can only be fulfilled if there is 
a single term in the sum. Thus a pure P-rep state, 
for which trp^ = 1, has to be a coherent state. 

6. Any density matrix can be decomposed as a sum 
of the totally mixed state Pq and a traceless hermi- 
tian operator p with trace norm one multiplied by 
a positive real parameter k, 

Pk = Po + «P- (8) 



Since C is convex, there is, for any given direction 
p, an extremal value Kg of n such that € C if 
< K < Kg and ^ C if K > Kg. The states 
p = pQ + KeP form the boundary dC of P-rep states. 
They belong to C provided we accept states p as P- 
rep if they can be approximated in the trace norm 
by a convex sum of coherent states, that is for all 
e > there exists a positive function P(a) such that 
/ daP (a) |a) (a|| I < e. With this extended 
definition the set of P-rep states becomes compact. 
In some directions the boundary dC may touch dJ\f, 
e.g. when p = \a) {a\ is a pure coherent state. 

7. dC is separated by a finite distance from the state 
Pq. In other words, all density operators in some 
finite neighborhood of Pq are P-rep. To show it 
let us choose P (a) containing only the mandatory 
components with K < 2j, 

P{a) = ^+P{(^), 

2j K 

Pia) ^ E PKQYKQia). (9) 

The Pkq are bounded since they are related to the 
coordinates pj^g of p by ([6]) and trp^ < 1. As the 
spherical harmonics are bounded on the sphere and 
|[9l) is a finite sum, there is an upper bound Pe to 
the non-trivial part P (a) when p and a are varied. 
Thus, all matrices pg -I- Kp with k < l/(47rPe) will 
be P-rep. 

III. P REP FOR SYSTEMS OF SMALL SPIN 

In the case of a spin-1/2 or a spin-1, it is possible to 
obtain a complete characterization of P-rep states. 

A. Spin-1/2 

We denote by cr —{ax,Uy,az) the vector formed by 
the Pauli matrices. Together with the identity matrix I2 
they form a basis of the space of 2 x 2 matrices. Any 
2x2 Hermitian matrix with unit trace can be written as 

p=i(l2+U.^), (10) 

and u is given by u = tr(pcr). The matrix p is non- 
negative if and only if |u| < 1. A physical density ma- 
trix p can thus be represented by a point inside the unit 
sphere (the Bloch sphere). Matrices corresponding to 
points on the unit sphere are pure states. Since for spin- 
1/2 any pure state is a coherent state, the convex hull of 
coherent states is the convex hull of pure states, which is 
the set of all density matrices. Thus all states are P-rep. 

It is straightforward to find an explicit decomposition 
in terms of angular momentum coherent states by simply 
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diagonalizing p, which leads to the sum of two projectors 
with two positive eigenvalues. Nevertheless, there is a 
large freedom in choosing the coherent states. Accord- 
ing to finding a P-representation for p amounts to 
finding positive real coefficients and projectors on co- 
herent states |ai)(aj| = i (I2 + n^'^.cr) with |n(*)| = 1 
such that p = J2i Since the (Ti form a basis of 

the 2x2 density matrices, this is equivalent to finding Xi 
and norm-1 vectors n^*^ such that 



(11) 



This can be trivially achieved e.g. by taking any pair of 
points on the Bloch sphere such that the line joining these 
two points contains the point representing u inside the 
sphere. 



B. Spin-1 

Let us now consider a spin-1 density matrix. We shall 
use the representation 



1 1 
P = -I3 ' 



1 \ JaJb + JbJa 
3— 2""" '2 ^ y ^b-^dabj ^ , 

a,b—x,y,z 

(12) 

where Ja are matrices of the angular momentum with 
j — 1. The Ja and the {JaJb + JbJa)/2, together with 
the identity matrix I3, form a basis of the vector space 
V of 3 X 3 hermitian matrices. Inverting relation lfT2|) we 
obtain 

Ua^tr{pJa), Wab ^T^r p {JaJb + JbJa) - Sab, (13) 

which shows that u gR'^ while is a 3 x 3 real symmetric 
tensor with trace 1. The projector on a coherent state 
In), written in the form (fT2l). reads 



is non-negative. 

Proof. First let us assume that the Eqs. lfT5|) do have a 
solution. Then Z can be written 



(17) 



and for any vector y g we have 

y*Zy = ^ A, (y.n(*)) ' - A.y.n^')^ > (18) 

since the weights A^ > sum to 1 and f{x) = is a 
convex function. Therefore Z is indeed non-negative for 
all P-rep operators p. 

Conversely, if Z > 0, then it is possible to exhibit 
a decomposition of p by finding an explicit solution to 
Eqs. dm). Let A be such that Z = AAK If we denote by 
t^') the eight column vectors (±1, ±1, ±1) obtained from 
all combinations of the ± signs, and define 




(19) 



then one can check that a solution to Eqs. ifTSj) is given 

by 

nW = u + T,At(*) (20) 



which proves that p is P-rep. 

The necessary and sufficient condition Z > in the 
case of spin-1 allows to characterize the boundary dC of 
P-rep states. Indeed, let us consider a one-parameter 
family of states as in |[8]) . If u and W are the vector and 
matrix corresponding to the expansion (fT2l) of the state 
Pq + p, then the vector and the matrix associated with 
Pk= Pq + i^P are given by 

u„ ~ KU (22) 



I \/ I ll , 1 Ti 1 ^ K \ JaJb + JbJa 

|n)(n| = -l3+-n.J+- [nanb--Sabj ^ ■ 



1 - K 



According to p is P-rep if and only if there exist 
Ai > with J2i = 1 and coherent states corresponding 
to vectors n^'^ G R"^ of length 1 such that 



(14) and thus the 3x3 matrix Z^ associated with reads 



1 - K 



I3 - K UU 



(23) 



(15) 



nl' = Wab, 



The value k — at which the scaled operator p^ ceases 
to be P-rep corresponds to the smallest n for which Z^ 
has a zero eigenvalue. Thus Ke is the smallest solution of 
the equation det Z^^ = 0, and the equation of dC in the 
vector space V is 



(with a,b running over x,y,z). It turns out that these 
equations admit a solution - and hence p is P-rep - if 
and only if the real symmetric 3x3 matrix Z with matrix 
elements 



1 



U = 1. 



(24) 



Zab = Wab - UaUb 



(16) 



This equation gives implicitly the value Kg for each di- 
rection p in the vector space V. As the examples of spin- 
1/2 and spin-1 show, the proportion of P-rep matrices 
among all density operators depends on j . 
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C. Necessary conditions for higher spins 



IV. NUMERICAL IMPLEMENTATION 



It is possible to derive more general necessary condi- 
tions for P-representability of spin-j states, as follows. 
Let us denote by Jt = t.J the spin operator in direction 
t. For a coherent state |n) corresponding to a vector n, 
the mean values of Jt and are given by 



(n| Jt|n) = i t.n 



(25) 
(26) 



Any P-rep state p can be written as p = 
^ ■ Ai|n('^)(n(*) I, which implies for the mean values of 
Jt and J^ in the state p 



( Jt) = j ^ Ai t.n 



(i) 



(Jt) = \ +3\3 



t.n 



(i) 



(27) 



(28) 



Convexity of /(x) 
to the inequality 



applied to the sums over i leads 



2j(Jt > - (2j - l)(Jt)' - / > Vt, |t| = 1, (29) 

with equality if and only if p is itself a coherent state. 
This is a necessary condition for P-rep, valid for any j. 
In the particular case of spin-1/2 this inequality becomes 
{Jt) ^1/4, which is obviously true for all states p and 



all directions t. In the case of spin-1 the inequality l[29|) 
can be rewritten as 

(2( J, Jb) - ( Ja) ( Jb) - (5afc) ta% > Vt = (i„ ty, t,), |t| 

a, 6 

(30) 

As can be seen from Eqs. (flSl) and (flGl) . this inequality 
exactly corresponds to the condition Z > derived in 
the previous section. 

For higher spins, one can similarly derive other neces- 
sary conditions. For instance for a P-rep state of spin- 
3/2, one has 



(j3).|5:A.(t.n(')) + 3^A.(t.n« 



(31) 



and a necessary condition imposed by the fact that 
I J2i '^i^il ^ J2t ^1^1 for any Xi € [-1, 1] reads 



A. Variational approach to P— representability 

Suppose we are given a density operator and want to 
establish whether it is P-representable. Let us use the 
multipole expansion The coefficients Pkq with < 
K < 2 j will be defined by Eq. ^ . Orthogonality of the 
spherical harmonics implies that the hypothetical P(a) > 
satisfies the integral equations 

J P(a)Y*KQ{a)da ^ Pkq, 0<K<2j, \Q\ < K, 



together with 



(33) 



J P{a)da = trp = 1. 



If we find any P(a) > satisfying these equations the 
state in question is P-representable. 

We can ask for more and try to find the representability 
boundary for all matrices of the form Pk, — Po + i^P ob- 
tained by scaling a given traceless normalized hermitian 
matrix p. To that end, we consider the set of matrices 
Pq/k + p, k > 0. These states all have the same traceless 
part p, thus they are represented by P-functions P(a) 
that satisfy Eqs. (|33l) with Pkq corresponding to p, but 
with / P(a)da — i. We look at the minimum of the 
functional F [P] = J P{a)da over these states. Suppose 
that the minimum is realized by some function Pe(a) and 
introduce Ke through 



min J P{a)da = J Pe{a)da = — . 



(34) 



_T^e corresponding density operator = Po + HeP is 
represented by the function KePe(a)- As we pointed out 
it means that all operators with < k < Kg are P- 
representable and that Pg belongs to the boundary dC. 



B. Concavity of I/ks 

The parameter Ke corresponding to the border of P- 
rep depends on the matrix p, such that Kg = Kg(p). Let us 
take two matrices, p^ and p^^ and calculate the respective 
Ke{p^), Ke{p^^). Consider now a convex combination 



cp' + c)p" , 0<c<l. 



Then 



Vt, 



< 



{J? 



(32) 



These necessary conditions can be considered as "non- 
classicality witnesses", as a state p is not in C if at least 
one of these conditions is not fulfilled. 



< 



«g(pW) - Ke{p') Keip")' 

i.e., 1/Kg is a concave function of p. The proof is based 
on Eq. (|34l) . Let Pg^ , P/^ be the functions minimizing 
J Pda under constraints corresponding to the operators 
p^ and p^^ respectively. Then the function P''^-' = cP/ + 



6 



(1 — c)P/^ will obey the constraints corresponding to the 
operator p^'^\ Therefore we must have 

— ^-TTT = min / P(a)da< f P^'^Ua) da 

= cj P/(a)da + (l-c) J P^' (a) da 
c _|_ -'^ ~ 

which impHes concavity of I/ks- Thus the knowledge 
of Ke for two density matrices gives a lower bound for 
a whole family of convex combinations of these density 
matrices. 



number is one less than predicted by Caratheodory's the- 
orem because the solution is a boundary, not an internal, 
point of the set of the density matrices P-representable 
by (|35|) . The minimum found numerically for a given 
n yields an upper bound on the exact value of l/^e 
(Eq. (|34l) ). i.e., the lower bound on the value of the scal- 
ing parameter k at the border of P-rep in p„ = Po + '^P- 
The linear programming approach was numerically 
tested and found efficient for moderate values of j. For 
a given p, the minimal value of diminished fast with 
the increase of n and was stable. On the other hand, 
the solution w changed erratically with the change of n. 
That was to be expected considering the freedom in the 
choice of P (a). 



C. Linear programming 

In order to numerically implement the variational ap- 
proach described here, let us choose the trial P-function 
in the form of a linear combination of (5-peaks 

n 

P (a) = ^ u;,^ (a - a,) (35) 

i=l 

where the points ai = {9i,4)i) are more or less uni- 
formly distributed on the unit sphere, and Wi > Q are 
non-negative variational parameters; the delta-functions 
are assumed to be normalized on the unit sphere, 
5{a — ai) = (5 (cos0 — cos^i) 5 ((/> — (/)j). Inserting this 
P(a) in (|33|) we come to the optimization problem: find 
w = {twi, . . . , Wn\ with all > 0, i = 1 . . .n, minimiz- 
ing the sum 

n 

P(w)=^w„ (36) 

1=1 

and subject to M = (2j + 1)^ — 1 linear constraints 

n 

J^YKQiai) ^ Pkq, 0<K<2j, \Q\ < K. 

This is a problem of linear programming [l^l- Its well- 
known theorem states that whatever the number of un- 
knowns n the minimum of F is realized on a solution 
containing no more than M non-zero components. This 



V. COMPOSITE SYSTEMS 

The definition of classicality can be extended to sys- 
tems of more than one particle in a natural way. In the 
present section we shall consider the case of two particles, 
but the formalism generalizes to an arbitrary number of 
particles. 

A. Classicality for two particles 

The P-representation of a density operator in the case 
of two spins jUand Jb, 

P"^ j d^aAd^aBP{aA,aB)\aA)\aB) {aA\{aB\ (37) 

with P > is possible for separable states only; con- 
sequently P-rep is a sufficient criterion of separability. 
The partially transposed matrices p"^-** and p^^ are de- 
fined in a fixed computational basis = ® \ j)B as 
pIj^m = Pk3,ii and p^^^^; = p.i^^^ . They are P-rep if and 
only if p is P-rep, and the corresponding P-functions 
P^-* and P^« are simply related to the P-function of 
p by P^-*(aA,aB) = P(aA,as), Q-a = {9 a, -'Pa), and 
correspondingly for P^^ . All previously considered equa- 
tions are reformulated for two spins in a straightforward 
manner; we shall list them without commenting. 

The representation of p in terms of products of spher- 
ical multipole operators reads 



I 

2jA Ka '^Jb Kb 

f = X! X! X! X! PKaQa.KbQb^KaQa^KbQb^ ^^^) 
Ka=0 Qa = -Ka Kb=0 Qb = ~Kb 

and we have the P-function expanded into products of spherical harmonics, 

oo Ka oo Kb 

P{(^)= X! X! X! X! PKAQAKBQBYKAQA{(^A)yKBQB{(^B)- 
Ka=0 Qa = -Ka Kb=0 Qb = -Kb 
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The relation between the coefficients of p and P is given by 

PKaQaKbQb = PkaQa.KbQb 



x47r- 



y/(2jA - KaWja + Ka + 1)!(2jb - Kb)1{2jb + Kb + 1)! ' 
and the density operator with a scaled non-trivial part by 

Po 



l(2j\4+l)x(2jfl+l) 

(2m + 1)(2jb + 1)- 



The following variational problem needs to be solved when the boundary of P-representability is to be found: 
minimize the functional 

F[P]= J SaASaBP{aA,aB) 
with P{aA,OLB) > satisfying the integral equations 

J (faA(faBP{aA, aB)Y^^Q^ (q:a)Y^^q^ (as) = PkaQa,KbQb > (39) 

where Ka,Kb run from to 2j excluding Ka — Kb ~ 0, and IQaI < Ka, \Qb\ < Kb- If the minimum of F is equal 
to 



Fe=minF= / (faAlfc(BPe{c(A,OiB) — — 
J Ke 



then the density operator lying on the boundary of P-representability will be p^^ . 

For the numerical implementation, the integrals are now taken over a product of two unit spheres of Alice and Bob. 
Let us choose the trial P-function as 

nA nB 

P{aA,aB)=^ ^w,^,^5{aA- af^)5{aB" af^) (40) 

iA = l iB=l 

where ua points and ub points af^ are uniformly scattered over the spheres of Alice and Bob, respectively, and 
WiAiB ^ cire UAnB variational parameters. We now solve the linear programming task: minimize 

riA riB 

"^lA iB 



^(w)=E E 

with WiA tB > satisfying M = {2ji + 1)^(2^2 + 1)^ — 1 Hnear constraints, 

nA nB 

E E ^KAQAi^tA)^KBQB(^fB) ^^AiB = PKaQa.KbQe 



iA=i iB = l 



Here Ka, Qa, Kb,Qb take all possible values excluding Ka = Kb = 0. Again, the optimal solution contains no more 



than M non-zero elements WiAiE 



B. Two spins 1/2 

Considering that the density operator of a single spin-1/2 is always P-rep it is easy to see that the density operator 
for a system of two spins is P- rep if and only if it is separable. Consequently, the necessary and sufficient condition 
of P-rep is given by the Peres-Horodecki theorem |25!, 2(1] . It means that the boundary of P-representability in the 
family — Po + ^P reached when either or its partial transpose p^* ceases to be non-negative. This was checked 
numerically in the linear programming approach: the minima l//te of the functional F [P] calculated with the matrix 
p and its partial transpose p"^'* in all cases coincided with each other and agreed with the scafing necessary to shift 
the smallest eigenvalue of either p or p^-^ to zero. The optimal P was obtained as a combination of M = 15 coherent 
states, some of them with very small weights. 
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FIG. 1: (Color online) Example of a set of classical states C for a bipartite system of two spins 1/2 and 1 parametrized by two 
parameters, p = Pq + i^^Pi + 1^292 with some traceless Pj^, Pj. Boundaries are shown of non-negativity of p (bold black line), 
non-negativity of its partial transpose p^^ (dashed line), and of P-representability of p, p^^ (inner red line). 



In this case the separabiHty and P-rep conditions do not coincide. Indeed consider for instance the pure product 
state (in \jm) notation) \^) = |ii)(g)|10). Then the mean value of the operator l2®Jf in the state IV') is (10| J^jlO) = 0, 
while using Eq. l|28|) one should have for a P-rep state (I2 ® J^) > 1/2. Thus, jV') is not P-rep. More generally, it is 
easy to show numerically that dC is well inside the separability boundary. An example is shown in FigdJ where we 
display the two boundaries for a density matrix of the form p = pg + nipi + K2P2 with two random but fixed traceless 
parts pi and P2- 



A simple necessary condition for P-rep can be formulated for the density operator p of the system of two particles 
A and B. Let Va be any non-negative operator in the Hilbert space of A and take the partial trace of pVa over the 
A— variables. Assuming that p is P-rep and using the coherent states \a') for the calculation of the trace we obtain 



C. Spins 1/2 and 1 



D. Classicality witness 




(41) 



(42) 



(43) 



where P {(3) = J daP{a,(3) (a|Vyi|Q!) is manifestly non-negative. Consequently, 

Pb = (Tr^ pVa) I TrpVA 



(44) 



can be considered as a density operator in the P— space which is P-representable by a function P (/?) / TrpV^i. 
Therefore p can be P-rep only if p^ is also P-rep (not vice versa) . The P-rep of p^ is easy to check using our result 
for J = 1. One can take, e.g., Va = 1a getting pg = Tr^ p. 
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VI. CONCLUSION 

The P-representable states are classical mixtures of 
projectors on angular momentum coherent states, i.e. of 
angular momentum states with minimal uncertainty. 
The P-rep states have many interesting properties. They 
can be seen as the "most classical" states, an "inner circle" 
within the linear space of density operators which forms 
a convex set C that contains the totally mixed state in its 
interior. In the case of two spins, C is a subset of the set 
of separable states. The study of the P-representation 
provides thus important information on the structure of 
space of density matrices. 

We have studied conditions for P-representability, and 
completely characterized the set of classical states for 
small spins: for a spin-1/2 all states are P-rep, and for 
a spin-1 we deduced a necessary and sufficient condition 
for P-rep. In the case of two spins-1/2, P-rep is equiva- 
lent to separability, but already for a spin-1/2 combined 
with a spin-1, there are states which are separable but 
not P-rep. In addition, we have shown that the question 
whether a given state is P-rep or not can be solved with 
a practical numerical method based on the linear pro- 



gramming algorithm for finding the border of P-rep. We 
have also formulated necessary conditions based on mea- 
surable observables for P-rep, which can be considered 
"non-classicality witnesses" for spin systems. 

Both analytical and computational methods have been 
used so far on very modest values of j (up to j ^ 2) ; for 
large j the numerical methods become forbiddingly slow. 
It would be important to investigate the limit of large 
j and provide thus a bridge to the case of continuous 
variables where the P-rep states were an object of intense 
studies for many years and proved to be of great physical 
importance. 
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